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Abstract: The program of understanding inverse-power law corrections to event shapes 
and energy flow observables in e^e~ annihilation to two jets and DIS (1+1) jets has been a 
significant success of QCD phenomenology over the last decade. The important extension 
of this program to similar observables in hadron collisions is not straightforward, being 
obscured by both conceptual and technical issues. In this paper we shed light on some of 
these issues by providing an estimate of power corrections to the inter-jet Et flow distribu- 
tion in hadron collisions using the techniques that were employed in the e^e~ annihilation 
and DIS cases. 
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1. Introduction 

The principle issue that limits the accuracy of theoretical predictions in QCD is the pres- 
ence of non-perturbative physics responsible for the confinement of quarks and gluons. The 
success of perturbative QCD, despite a lack of quantitative understanding of the confine- 
ment process, is a major achievement that is based on identifying infrared and collinear-safe 
observables [1] which are as insensitive as possible to non-perturbative effects such as hadro- 
nisation. In such instances the role of non-perturbative effects is reduced to the level of 
corrections to the perturbative estimates which take the form of an inverse-power law in 
the hard scale Q, where p depends on the observable. For some observables, such as 

the total cross-section in e~^e~ annihilation to jets, where p = 4, the corrections in question 
are insignificant and can safely be ignored in comparisons to data. For other observables 
such as event-shape variables (see ref. [2] for a review) it was noted however that the power 
corrections scale as 1/Q and obscure the perturbative analysis significantly. 

Over the past decade, theoretical efforts essentially based on renormalons (see ref. 
[3] and references therein) have given a clearer picture of the origin and role of power 
corrections. Within the renormalon model these corrections are shown to be related to a 
factorial divergence of the perturbative QCD expansion at high orders and the consequent 
error in truncating what is in fact an asymptotic series [3] . This observation allows one to 
estimate power corrections from a lowest-order Feynman graph, modified to incorporate 
the relevant subset of higher-order terms (renormalon bubble insertions). 

From a phenomenological viewpoint, perhaps the most widely used formulation of the 
renormalon model is one that uses a dispersive representation of the running coupling 
[4], thereby introducing a dispersive variable m that plays the role of a fake gluon mass 



- 1 - 



[5]. This scale is a natural trigger for power corrections arising in the infrared regime 
for observables that arc otherwise dominated by a hard scale Q. The most appealing 
feature of the dispersive approach is the hypothesis of a universal infrared-finite coupling 
Ug. The overall size of the power corrections in this approach is related to moments of 
this universal coupling: dm^ lm?as{rii?){m? [Q'^Y . For most event-shape variables (for 
example) p = 1/2 [2], and the relevant coupling moment or equivalently the related quantity 
ao has been extracted from data for different event shapes in both e'^e~ annihilation 
and DIS [6,7]. The values of ao thus obtained have generally confirmed the universality 
hypothesis to within the expected uncertainties due to missing higher-order corrections [2] . 
These studies have not only enabled successful studies of several observables, but also lent 
credence to the notion that the QCD coupling may be a meaningful (finite and universal) 
concept all the way down to the smallest energy scales, and have thus renewed hope of a 
better understanding of the confinement domain from first principles of QCD. 

The success of renormalon-inspired studies has thus far been limited to observables 
which involve just two hard partons at Born level such as event shapes in e+e~ — two 
jets and DIS (1 + 1) jets. Some of the most interesting QCD observables however do not 
fall into the above class. Examples include three-jet event shapes in e^e~ annihilation and 
DIS [8-11], dijct event shapes in hadron collisions [12] and single-jet inclusive cross-sections 
at hadron colliders [13]. These observables contain more than two hard partons at Born 
level and additionally pertain to processes that have gluons at Born level where one may 
question the extension of the dispersive approach, which (strictly speaking) was formulated 
for observables involving only quarks at the Born level [4] . The program of understanding 
power corrections has not been put to test in such situations which will be of importance 
at the LHC, for instance, although some progress is being made on the phenomenological 
side for three-jet event shapes [11]. In fact for the case of observables involving four hard 
partons at the Born level there are as yet no theoretical predictions for power corrections. 

Applying the renormalon model, in many such cases, reduces (as for e^e" annihilation 
to jets or DIS event shapes) to an analysis of soft gluon radiation with transverse momenta 
kt ~ AqcD) which are associated with a universal infrared- finite coupling. Given the rich 
structure (colour and geometry-dependence) of soft gluon radiation for processes like dijet 
production in hadron collisions, it is indeed enticing to see how a perturbative structure 
may influence predictions for power corrections as predicted by the renormalon model. 
Signs of perturbatively calculable colour structure in the non-perturbative power-behaved 
component would strongly suggest that the renormalon-inspired picture of these corrections 
is an extension of soft gluon radiation, with a modified but universal coupling, thereby 
conclusively establishing the model. 

In this paper we provide a calculation of the power correction to the transverse energy 
{Et) flow away from jets accompanying hard dijet production in hadron collisions. A variant 
of this observable (the "pedestal height") was suggested several years ago as a means of 
separating perturbative bremsstrahlung from the contribution of the soft underlying event 
in hadron collisions [14]. To date however a satisfactory understanding of the away-from- 
jet energy flow has proved elusive. In the region of small Et the distribution da/dEt 
contains large logarithms of the form d^\n^~^{Pt/ Et)/ Et, where Pt is the hard scale (jet 
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transverse momenta) of the problem. These logarithms can only be resummed in the large 
Nc limit [15-17] which limits the accuracy of perturbative estimates that can be made 
in this case alongside the lack of any estimates of the NLL contributions. One may also 
expect power corrections of the form 1/Pt to the resummed distribution as for the case of 
event-shape distributions in e+e~ annihilation and DIS. Once again there is no estimate 
for these power corrections and it is this that we aim to provide here. Combining the 
power corrections we compute with resummed Et flow distributions gives a more complete 
theoretical account which should facilitate comparisons to experiment. 

While at hadron colliders the ever-present soft underlying event may obstruct clean 
studies of power-behaved corrections arising from the bremsstrahlung component of the 
Ef flow, our calculations here are also easily adapted to the case of rapidity gaps in dijet 
photoproduction at HERA, and are in principle more readily tested in that environment. 
Moreover there are theoretical ideas concerning obscrvablcs such as inclusive jet cross- 
sections at hadron colliders, where it may be possible to disentangle the underlying event 
from the non-perturbative physics in the bremsstrahlung component, due to a singular 1/R 
{R being a jet radius parameter) behaviour of the latter [13]. A full estimate of this piece 
however involves a calculation similar to the one we introdTice here and thus the present 
calculations can be used as a guide in moving towards better estimates of the inclusive jet 
cross-sections as well. Another related class of observables, to which the techniques we use 
here are directly applicable, is the important case of event shapes at hadron colliders [12], 
for which resummed perturbative estimates already exist and identical considerations to 
those of this paper will be required when dealing with the issue of power corrections. While 
other problems such as lack of knowledge about next-to-leading logarithms also need to 
be addressed, the calculation of the power correction is ultimately an important ingredient 
which, as we stressed above, also serves as a case study for hadron-hadron observables more 
generally. 

This paper is organised as follows. In the following section we define the observable 
more precisely and review the perturbative result in Mellin space conjugate to Et, which 
involves colour matrices in the resummed anomalous dimensions [18-21]. We then compute 
the power correction to each of the matrix elements of the anomalous dimension by consider- 
ing the appropriate combination of dipoles involved in that matrix element. The calculation 
of the power correction in each dipole term is performed using the appropriate scale of the 
running coupling (the invariant transverse momentum of the dipole, k± [8-10,22-24]). The 
final step is to take the inverse Mellin transform of the result including non-perturbative 
corrections, for which it proves convenient to diagonalise the power-corrected anomalous 
dimension matrices. We find that the power correction to the Et distribution is not a simple 
shift of the resummed distribution by a fixed amount proportional to 1/Pt as is the case 
for event shapes in e+e~ — > two jets and DIS (to a good approximation). 

2. Resummed perturbative result 

Here we outline the resummed result for the Et flow distribution accompanying hard dijet 
production in hadronic collisions, which was originally computed in ref. [25], specialising 
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for simplicity to the case of a slice in rapidity of width Ai], which we take to be centred at 
7/ = (the rapidity rj and transverse energy Ef are both defined with respect to the beam 
direction). We also work, again purely for the sake of convenience, with fixed jet transverse 
momenta Pt and rapidities but can easily generalise to other differential distributions of 
the hard dijet system. The observable we compute is then the integrated distribution: 



where o"o is the Born cross-section for dijet production (at fixed jet transverse momenta and 
rapidities) in hadron collisions, O denotes the rapidity slice and Et is obtained by summing 
over all objects (partons/jets) in the gap. 

Since we are examining the accompanying Ef flow distribution, we are dealing typ- 
ically with Et <C Pt and hence there arises the need to resum large logarithms of the 
form a'^ In"^ (Pt / Et) , with m < n. The resummed prediction for Ti{Pt,Et) up to single- 
logarithmic accuracy (resumming all the terms a" ln"(Pt/£'t) in the logarithm of the in- 
tegrated distribution), can be thought of as comprising two distinct pieces with different 
dynamical origins: "global" and "non-global"^ [15,16]. 

The first piece (the so-called "global" component) is a result of considering multiple 
soft emissions, both real and virtual, which are attached just to the primary or Born hard 
partons. Due to infrarcd-safcty of the observable real and virtual emissions cancel below 
the scale Et, while real emissions above this scale arc vetoed. Thus the resummed result 
for this piece is just the summation of virtual graphs above the scale Et, attached to the 
primary hard partons [18-21]. In actual fact the factorisation of real emissions and the 
consequent cancellation with virtual ones takes place in Mellin space conjugate to Et. This 
complication can be ignored for the leading single-logarithmic terms, but to analyse the 
impact of power corrections on the resummed distribution we need to compute the result in 
Mellin space and then invert the transform to Et space. To be more precise the perturbative 
resummed result (considering just the above described global term for now) reads (see for 
instance ref. [17]): 



where v is the Mellin variable conjugate to Et and the integration contour is taken, in 
the usual manner, parallel to the imaginary axis and to the right of all singularities of the 
integrand. We have: 



In the above V is essentially an "anomalous dimension matrix" and H and S are the 
"hard" and "soft" matrices. The matrix elements Hij represent the product of the Born 
amplitude in colour channel i and its complex conjugate in colour channel j, and the 
matrix S represents the normalisation arising from the colour algebra. The squared matrix 

^ We ignore, for now, a possible additional complication afflicting non-global observables in hadronic dijet 
production - super-leading logarithms [26] that may appear at fourth-loop level. The status of these is as 
yet unclear. 
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element for the Born scattering in this notation is just = Tr(HS). We shah return to 
the detailed structure of T in the next section but it is important to point out here that all 
the matrices above depend on the partonic 2^2 sub-process considered, e.g. they differ 
for say qq — > qq and gg gg sub-processes. Their form also depends on one's choice of 
colour basis. For examples the reader is pointed to the original references [18-21]. We also 
point out that since we deal exclusively with soft radiation and differential distributions in 
the jet rapidities and transverse momenta, the parton distribution functions simply factor 
from the resummed result and cancel against those in the normalisation factor ctq. 

We now turn to the second piece of resummed distribution: the "non-global" contribu- 
tion. The above result, resumming essentially virtual corrections above the veto scale Et 
which dress the hard scattering, is not the complete description at single-logarithmic accu- 
racy [15,16]. An additional "non-global" piece S{Pt,Et) arises (starting at the two-gluon 
emission level, 0(0^)) and is also single- logarithmic. The dynamical origin of this piece 
is multiple soft energy-ordered radiation with an arbitrary complex geometry ( "hedgehog" 
configurations of soft gluons, as opposed to emissions strongly ordered in angle). It has 
thus far been possible to treat this term only in the large Nc approximation, which limits 
the accuracy of perturbative calculations in the present instance. 

However it has been pointed out and clarified in a series of papers [27-30] that the role 
of the non-global component can be significantly reduced by defining the observable in terms 
of soft jets rather than individual hadrons in the gap. This can be achieved by running 
a clustering algorithm on the final states such that all objects are included in jets. The 
clustering procedure (with a large cluster radius R = 1) was shown to virtually eliminate 
the non global component while giving rise to additional global terms [29] that were at most 
modest corrections to the pure virtual dressing as represented here by eq. ( ^ ) . Moreover 
the power corrections associated to the non-global component of the result would start at a 
higher-order in (albeit potentially accompanied by logarithms of Pt/ Et) which we ignore. 
This was also the procedure employed for the case of non-global DIS event shapes [2] , where 
power corrections were computed in the exponentiated single-gluon piece of the resummed 
distribution, which was a phenomenological success. For all these reasons we shall choose 
to concentrate for the rest of this paper on the global component (eq. (2.3)), ignoring the 
non- global component. 

In the following section we explicate the structure of T in terms of the various hard 
colour dipoles from which one considers soft gluons to be emitted according to the usual 
antenna pattern. Non-perturbative power corrections are then computed on a dipole- 
by-dipole basis, adapting the procedure for a qq dipole developed for the case of e^e~ 
annihilation to two jets. Having obtained the power corrections to V[v) we then invert the 
Mellin transform to examine the result for the Et distribution. 



3. The anomalous dimension and power corrections 

We first write down the structure of the resummed anomalous dimension matrix T{v) 
and then note that it contains an integral over the running coupling which is formally 
divergent. Making the ansatz of a universal infrared-finite coupling cures this divergence 
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and introduces calculable power corrections to the perturbative anomalous dimensions. In 
what follows and for the rest of this paper we specialise to the case of the sub-process 
qq — > qq since identical considerations are involved for all other sub-processes. Full results 
combining all sub-processes will be made available in forthcoming work. 

For the sub-process q{pi,ri) + q{p2,r2) 9(^3,^3) + q{p4,r4), where pi and are 
respectively four-momenta and colour indices, we can choose to work in the f-channel 
singlet-octet colour basis: 



For this basis we have: 

^=(Z^ nIi], (3-2) 



and the anomalous dimension matrix T is: 

\S-U CFS-^^{T-2U + S)j ' 

In the above S, T and U arc combinations of dipole contributions with each contribution 
given by the corresponding dipole antenna. Thus one has^: 

S =wi2+ W34., (3.4) 
T =wis + W24, (3.5) 
U = ■f2;23 + 'i^i4, (3.6) 

where each dipole contribution Wij reads: 

f d?k 2 {Pi-Pj) lu X .Q'7^ 

In the above result we have integrated the soft gluon emission probability, given by the 
dipole antenna pattern, over the gluon phase-space (with k and w being respectively the 
three-momentum and energy of the gluon, as usual) with a "source" function u{kt., and 
g1 = Airag. The source is a result of factorising the real soft emission phase-space in Mellin 
space (see for instance ref. [17]) and accounting additionally for virtual corrections. It 
reads: 

u{h,u) = (l-e-'''^), (3.8) 

if the emission is in and zero elsewhere. The source thus represents the impact of 
real-virtual contributions which completely cancel, to our accuracy, for emissions outside 
O. Now introducing the variables r] and 0, respectively the rapidity and azimuth of the 
emission with respect to the beam direction, one can write: 



^Here our use of the 5", T and U labels differs slightly from that of other references, e.g. [18-21]. 
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with fij{ri,4>) being the functional dependence on rapidity and azimuth that arises from 
the dipole antenna patterns and which we shah use below (see appendix). We note that, 
to our accuracy, the correct argument of the running coupling for emission from a dipole 
( [8-10,22-24]) is k'j_-j = 2{pi.k){pj .k)/{pi.pj), which is the transverse momentum of the 
gluon k with respect to the dipole axis in the dipole rest frame. This must be distinguished 
from kt, the transverse momentum of the gluon with respect to the beam direction, which 
is the quantity that directly enters the observable definition. In fact we have, in terms of 
the functions fij introduced above, k±^ij = kt\j2jjij. 

3.1 Power corrections dipole-by-dipole 

Now we proceed to an extraction of the leading power-behaved contribution. In order 
to do this we first note that the integral over kt in eq. ( |3.9| ), which can be rewritten as 
one over the related variable is divergent if one uses the usual perturbative definition 
of a^, due to the divergence of the perturbative running coupling at k^ = Aqcd- In 
order to isolate and cure this pathological behaviour we assume an infrared-finite coupling 
{as) and change variable of integration from kt to k± in eq. ( p. 9] ). We then follow the 
method of ref. [31] to write as{k\) = as,pT(^D + '^Os,np(^D) where PT and NP stand 
for perturbative and non-perturbative respectively. In doing so we have assumed that the 
actual coupling a<j is in fact finite even at arbitrarily small k± , and can be split into the usual 
perturbative component a^^pT and a modification Ja^^NP which is due to non-perturbative 
effects. Both the perturbative and non-perturbative components separately diverge, but 
the divergences cancel in their sum due to the assumed finiteness of the physical coupling 
Ug. Moreover, since we do not wish to modify the perturbative results at large scales, 
the non-perturbative physics as represented by the modification <5as,NP niust vanish above 
some infrared "matching" scale /i/. Effectively the addition of the 5a<j^NP term represents 
removal of the badly-behaved perturbative contribution below fij and its replacement with 
the well-behaved integral over the infrared- finite physical coupling a^. 
Thus for the observable itself one has from dipole (ij): 



The integral involving the perturbative coupling represents the usual perturbative contri- 
bution from dipole (ij). The leading logarithmic perturbative contribution arises from the 
region where one can make the approximation: 



The perturbative results are reported at length in ref. [25]. In what follows we shall consider 
in more detail the non-perturbative contribution from the integral involving (Ja^^NP- 

In order to evaluate the non-perturbative contribution we first consider that the leading 
such term arises from the region /i/ <C l/i^, which translates to a requirement on Et to be 




(3.10) 
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above a few GeV. In this region one can expand the exponential in an exactly analogous 
way as for event-shape distributions [31]. The leading term is given by the first term in the 
expansion: 1 — exp(—ukt) ~ i^kt, and this corresponds to a linear 1/Pt power correction. 
We ignore quadratic and higher power corrections that would scale as 1/Pf and beyond, 
once again following the case of event-shape variables. We also note that in the shape 
function approach [32-34], where one may study non-perturbative effects even into the 
region Et ~ Aqcd, higher powers of need also to be retained. Working with just the 
leading term gives us the non-perturbative correction from the (ij) dipole which can then 
be written as: 

Wij = wf^ + uVaj. (3.12) 

Here the non-perturbative quantity V is the first moment of the coupling modification 
(5qs,np: 



which also enters 1/Q {Q being the hard scale) power corrections to event shapes and can 
be related to the parameter ao, extracted from fits to event shape data, as: 

V = ^ (ao(w) - asiP?) - ^ (in ^ + f + l) «'(^') + ^(«')) ' (3-14) 

where ao(/^/) = l/w Jq' dk± as{k'j_), /3o is the first coefficient of the QCD beta function 
and: 

The coefficients Cij represent the integral over directions: 



where fij arises from the dipole antenna pattern as indicated in eq. ( |3.9D , and a further 
factor proportional to y^Jij comes from rewriting kt in terms of k± as we stated before. 
The explicit form of the fij functions is reported in the appendix. Performing the integrals 
over rj and (j) in eq. ( |3.16| ) yields the coefficients Cij that correspond to the non-perturbative 
contribution to Wij in eq. (l3.12| ) , which we do not explicitly display for economy of presen- 
tation. 

Having computed the power corrections proportional to v for each dipole, we can 



include these corrections to the anomalous dimension matrix in eq. (3^), which can then 
be written as: 

r = T{u)TpT + i^VTf^P, (3.17) 

where the non-perturbative contribution Fnp is built up by combining the dipole contribu- 
tions Cij as in the perturbative case. In the case of the perturbative term Tpx we explicitly 
extracted the integral over the transverse momentum of the coupling: 

dk, a.,PT(t?) 
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which arises by making the substitution ( 3.11 ) in the first term of ( 3.10| ). Then the matrix 
FpT is the usual perturbative anomalous dimension containing integrals over gluon direc- 
tions inside the region^ Q. In the following section we shall consider the evaluation of the 
inverse Mellin transform to take our results from i/ space to Et space. 



4. Power corrections in the Et cross-section 

After accommodating the leading power corrections (those expected to give rise to l/Pt 
effects) eq. ( |2.2D assumes the explicit form: 



2711^ 



X TV 



Hexp (-tT\,^ - zypr};jp) S exp (-rFpT - vVT^^)] /Sq. (4.1 



In order to invert the Mellin transform (perform the v integral above) it is simplest to 
diagonalise the matrix tFpt + vVV^-p ■ In the basis in which the matrix V is diagonal the 
matrices H and S become H = R^^HR"^^^ and S = R^SR, where R is a matrix which 
contains the eigenvectors of T as column entries (see also [35]). After diagonalisation we 
can write the result for S(Pf, Et) in terms of components as: 

where are the eigenvalues of the matrix T. For the case of the qq — > qq sub-process, 
which we use as an example, F, H and S are 2x2 matrices and the above result can be 
written explicitly in terms of the elements of the various matrices as: 



nPuEt) = j 



dv ,,rp 1 



27riz^ So 



where the above result contains both perturbative and non-perturbative contributions. To 
separate these we note that the eigenvalues can be expanded so as to retain only the 
first-order in correction to the perturbative value, which depends logarithmically on i^: 

Ai = r(i^)Ar + z.PAfP + 0(z.2). (4.4) 

We emphasise here that while Xf^ are simply the eigenvalues of FpT, A^^ are not the 
eigenvalues of Fnp- Instead they are coefficients of the 0{i') component of the expansion 
of the eigenvalues of rFpT -|- z^T^Fnp and they depend on the components of both Fnp and 

FpT- 



These integrals are similar to those which yield the dj except that the functions fij are involved rather 
than /f/VV2. 
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The matrices H and S also differ from their pure perturbative forms by corrections 
which depend on v. Let us introduce the notation Dij = HijSji (summation not imphed). 
We expand the elements Dij to first order in f and write the result as: 

^(^.■^.) = E / ^eK-.-K'.r'i) (a- ^ ,-.M(Ar-«.",/^„, 

(4.5) 

where the sum runs over the components of the matrix Dij and we ignored ©(z^^) terms in 
the expansion of Dijiv), as they would contribute to corrections that scale as which 
are beyond our control with the present method. We shall see that the terms we ignore 
are expected to be numerically of no significance. 

We first write down the pure perturbative result Ti^'^{Pt, Et) obtained by ignoring all 
non-perturbative (NP) components. To perform the i' integral we use the standard method 
of expanding t{u) about the saddle point at u = l/Et, 

.M = r(W)+M.P,)(^)^^^^^_ + .... (4.0) 

Performing the integral over and noting that the logarithmic derivatives of t(z^) are 
related to subleading terms of relative order to the leading single logarithms, we arrive 
at the leading-logarithmic resummed result: 

SPT(P„ E,) ^ Yl D^e-('^^^'*r>^^^^/^o, (4.7) 

where the effect of the integration amounts to merely replacing by 1 /Et . Now redoing 
the integral including the non-perturbative correction to the eigenvalues (still ignoring 
the NP corrections to Dij) we find by examining eq. (4^) that the impact of the non- 
perturbative term amounts simply to a shift of the perturbative result in each of the terms 
in the sum in eq. (|4.7| ): 

i?i-i?i-P(AfP + AfP). (4.8) 

Looking at the distribution in Et, with Et being measured in units of the hard scale Pt, 
amounts to a 1/Pt non-perturbative shift in each term in the sum above, as is the case 
for two-jet event-shape variables [31,33]. However in contrast to the case of event shapes 
it should be clear that the overall impact of the power correction is not simply a shift of 
the perturbative distribution by a fixed amount since each term in the sum on the right 
hand side of eq. ( [4.7[ ) receives its own characteristic shift depending on the sum of the 
eigenvalues A^^ -|- A*'"^^ entering the term in question. 

We have still not accounted for the non-perturbative contribution to the colour basis 
as contained in the D^j^ terms. To evaluate these one performs the contour integral in 
question which yields a power correction of the form V/Et, which is related to the fact that 
these pieces are proportional to v. Computing the full result for this piece using precisely 
the same expansion about v = l/Et as before and discarding perturbative subleading terms 
beyond single-logarithmic accuracy, we arrive at our final result for Ti{Pt,Et): 



E(p„ii;,) = 5;e-^SM^- 



PT , nNP ^ 1 ^ / _ Pt 



/So, (4.9) 
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Figure 1: Power corrections to tlie energy flow distribution. PT stands for the pure perturbative 
result S^'^(Pt,£'t) as presented in eq. ( [4.7| ), PT+NP stands for the result including the non- 
perturbative correction and ignoring the terms Df,^, while PT+NP2 stands for the result presented 



in eq. (4.9) 



where A^P = viXf^+Xf^) and Af^ = \fT + xf^, and the function dj is approximately 
a constant of order as{Pt), varying very slowly with Et over the range of Et we consider 
here. It is a function of the logarithmic derivative of the single-log resummed perturbative 
result and hence scales as a,, 



G 



d 



Afj^ T{Et) +Hr{Et)) ^ as{Pt). 



(4.10) 



dln{Pt/Et) 

We also note the presence of l/r(£'t) accompanying the 1/Et dependence above, which is 



a reflection of the fact that the correction terms (involving Df-^) go as v/t. 



5. Results and conclusions 

In this section, for completeness, we illustrate the impact of non-perturbative power cor- 
rections on the energy flow distribution we discussed above. 

In fig. H we show the result for Ti{Pt,Et). We present the pure perturbative result 
(eq. (|4.7|)), the result including non-perturbative corrections without the DfF component 



(PT+NP) and the result presented in eq. (|4j) (PT+NP2). 

The results above were obtained for the illustrative value of pt = 80 GeV and fixing 
jet rapidities at -2.5 and 0.9 units respectively. We have performed the integration over 



the functions f^-"^ for the non-perturbative component (given in eq. ( p. 16 )) numerically. 



We also assumed that the rapidity gap has width A?? = 1. 

We notice that the effect of the term Df^ is in fact very small and can be safely 
ignored which is as an indication that neglected higher orders in the expansion of Dij 
are even more suppressed. We obtain as expected only a few percent correction to the 
resummed perturbative result even at fairly low Et, which given the remaining theoretical 
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uncertainties (notably having only the large Nc control over even leading logarithms) is 
not significant in itself. However our main motivation as stressed in the introduction, has 
not been the immediate phenomenology of the Ef flow observable, but rather a study of 
how power corrections may be included in resummed predictions for general observables 
in hadronic dijet production, since the soft function we treated here generically appears in 
resummation of global hadronic dijet observables. 

In conclusion, we have carried out for the first time a calculation of power corrections 
to transverse energy flow associated with dijet production in hadron-hadron collisions. We 
illustrated the computation for the energy flow distribution using the process qq qq. 
The generalisation to other channels is straightforward and only requires the numerical 
computation of the diagonaliscd anomalous dimensions and the corresponding hard and 
soft matrices. We found that the result does not correspond to the usual shift found 
in studies of two-jet event shapes and energy flows. The reason for this is the non-trivial 
colour algebra involved in the case of hadron collisions. The techniques we used here should 
enable better estimates of power corrections for observables which have a similar nature 
to the one we introduced here, such as the global inclusive jet cross-section we mentioned 
earlier. 

We re-emphasise that the simple calculation of this paper is just a first step in the full 
quantitative estimate of power corrections, pending the inclusion of two-loop effects. In 
the simpler cases of e+e^ 2 jets and DIS (H-l) jet processes, two-loop effects merely 
rescaled the simple one gluon calculation by a universal (Milan) factor [36-39] , which arises 
when considering a two-loop analysis for the argument of the coupling. We shall leave the 
inclusion of such effects to forthcoming work. 
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A. The functions 

We present here expressions for the functions fij. We define % and 774 to be the rapidities 
of the outgoing hard legs. We specify the kinematics of the particles as follows: 



Pi= (1,0,0,-1), (A.l) 

P2 = ^2^^/2(1,0,0,1), (A.2) 

P3 = -Pt (cosh ?73, 1,0, sinh%) , (A.3) 

P4 = Pt (cosh 774, —1,0, sinh 774) , (A. 4) 

k = kt (cosh rj, cos 4>, sin (f>, sinh r]) , (A. 5) 



where s is the hadronic centre-of-mass energy squared, related to ,s (the partonic centre- 
of-mass energy squared) by s = X1X2S, with xi and X2 being the momentum fractions of 
the incoming protons, carried by the incoming partons "1" and "2" respectively. 
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The functions fij{'q,<p) = pi.pj /{pi.kpj.k), are given by: 
/i2 = 2, 



/l3 
/l4 
/23 
/24 
/34 



cosh(?73 — rj) — 


cos </) 






COSh(T/4 — ??) + 


COS ^ 






cosh(?73 — 7]) — 


COS 4> 







cosh(?74 — r?) + cos <j) ' 

cosh(?73 - 174) + 1 
(cosh(r73 — rj) — cos (/))(cosh(r/4 — r/) + cos < 



(A.6) 
(A.7) 

(A.8) 

(A.9) 

(A.IO) 

(A.ll) 
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